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ABSTRACT 


The  greatest  of  constant  false  alarm  rate  processor  (GO  CFAR)  is  a  useful 
architecture  for  adaptively  setting  a  radar  detection  threshold  in  the  presence  of  clutter 
edges.  The  GO  CFAR  input  is  often  the  envelope  detected  in-phase  (1)  and  quadrature 
(Q)  channels  of  the  baseband  signal  (xe  =  sqrt  (I2  +  Q2)).  This  envelope  detection  can  also 
be  approximated  using  x=a|I|+b'Q|  which  requires  less  complex  hardware  (a  and  b  are 
simple  multiplying  coefficients).  The  envelope  GO  CFAR  processor  and  several 
envelope  approximation  GO  CFAR  processors  are  compared  in  terms  of  the  probability 
of  detection  (PD)  performance.  Closed-form  expressions  which  describe  the  PD 
performance  are  given  and  their  accuracy  evaluated. 
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I.  INTRODUCTION 


A.  GO  CFAR  PROCESSOR 

The  greatest  of  constant  false  alarm  rate  (GO  CFAR)  processor  is  a  useful 
architecture  for  adaptively  setting  a  radar  detection  threshold  in  the  presence  of  clutter 
edges.  The  inputs  to  the  processor  often  involve  the  envelope  detected  in-phase  (I)  and 
quadrature  (Q)  channels  of  the  baseband  signal  xe  =  JT2  +  Q2  .  This  envelope  detection 
can  be  approximated  using  x  =  amax{|/|,|£?|}  +  6min{|/|,|£?|}  which  requires  less 
complex  hardware  and  where  a  and  b  are  simple  multiplying  coefficients  (Pace,  1994 
and  Hache,  1994).  This  approximation  can  be  simplified  further  to  x  =  a\l\  +b\Q\  which 
reduces  hardware  requirements  even  more  That  is.,  the  two  channels  do  not  have  to  be 
compared  to  determine  the  minimum  and  the  maximum.  Ironically,  this  reduction  in 
complexity  has  increased  the  difficulty  of  obtaining  closed  form  analytical  solutions  to 
the  GO  CFAR  performance  in  terms  of  probability  of  false  alarm  (PFA)  and  probability 
of  detection  (PD).  Recent  studies  have  used  numerical  and  Monte  Carlo  methods  to 
detail  the  performance  of  this  processor  in  terms  of  the  PFA  as  a  function  of  the 
threshold  multiplier.  Curve  fit  solutions  have  also  provided  closed  form  expressions  for 
the  PFA  (Pace,  1994)  Closed  form  expressions  for  the  detection  performance  however, 
have  not  yet  been  examined. 

B.  PRINCIPLE  CONTRIBUTIONS 

This  paper  extends  the  results  from  previous  studies  and  develops  closed  form 
expressions  for  the  PD  performance  as  a  function  of  the  signal-to-noise  ratio  (SNR). 
These  results  provide  system  designers  with  an  accurate  estimate  of  system  performance 
for  several  detector  approximations  without  requiring  numerical  analysis  or  Monte  Carlo 
simulations.  Since  no  closed  form  expressions  for  the  PD  exist,  Monte  Carlo  simulations 
are  first  used  to  obtain  a  full  set  of  the  PD  curves  using  a  PFA  of  KT4  for  the  number  of 
reference  cells  n=  1,2, 4,8. 16.32  and  for  various  values  of  multiplying  coefficients  a  and  b. 
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Various  curve-fit  techniques  are  then  used  to  obtain  closed  form  expressions  for 
the  PD  vs.  SNR  curves  Several  formulations  are  investigated.  The  most  promising 
expression  involves  the  use  of  the  cif  function  which  produces  a  detection  performance 
within  ±0.015  of  the  actual  PD. 

Next,  the  coefficients  for  each  envelope  approximation  PD  performance  are 
plotted  as  a  function  of  the  number  of  cells  n.  This  allows  the  extraction  of  the 
appropriate  coefficients  for  any  number  of  reference  cells  n  between  0  and  32.  This  in 
turn  provides  a  quick  and  relatively  accurate  expression  for  the  PD  as  a  function  of  the 
SNR  for  the  various  detection  approximations  a  and  b.. 

C.  THESIS  OUTLINE 

Chapter  II  starts  with  a  review  of  basic  signal  detection  techniques,  an 
introduction  to  the  GO  CFAR  processor,  and  the  concepts  of  noise  and  probability.  It 
continues  with  a  description  of  envelope  detection  as  well  as  several  approximations  to 
the  envelope  detector  and  concludes  with  an  analysis  of  the  GO  CFAR  processor 
performance  in  terms  of  the  PFA  and  the  PD.  Chapter  III  provides  a  discussion  of  the 
Monte  Carlo  techniques  and  the  simulation  results.  Chapters  IV  presents  the  derivation  of 
closed  form  expressions  for  the  detection  performance  and  provides  tables  of  the 
corresponding  coefficients.  Chapter  V  presents  the  plots  of  the  closed  form  expression 
coefficients  as  a  function  of  the  number  of  cells  n.  The  appendix  contains  the  code 
developed  for  this  research. 
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II.  CFAR  PROCESSOR 


A.  BASIC  SIGNAL  DETECTION  PRINCIPLES 

The  detection  process  for  a  typical  radar  system  involves  an  inherently 
straightforward  comparison  between  the  received  signal  and  a  reference  value.  In  a 
simple  detector,  the  reference  value  is  a  fixed  threshold  which  remains  unchanged 
through  the  detection  process  If  the  received  signal  is  above  the  threshold,  a  target  is 
assumed  to  be  present.,  otherwise  the  signal  is  ignored.  Not  surprisingly,  this  is  referred  to 
as  threshold  detection.  Figure  1  illustrates  this  concept. 


Typical  Envelope  of  a  Radar  Receiver  Output  as  a  Function  of  Time.  (After  Skolnik,  1980) 

A  low  threshold  increases  the  chance  that  noise  alone  will  rise  above  the  threshold  and 
generate  false  alarms.  On  the  other  hand,  a  high  threshold  will  result  in  missed  targets. 
Consequently  the  selection  of  a  proper  threshold  is  a  compromise  between  reducing  the 
PFA  and  increasing  the  PD. 

To  make  things  more  difficult,  subtle  changes  in  the  threshold  can  drastically 
change  the  overall  system  PFA.  Although  these  changes  can  be  manually  compensated 
by  trained  operators  (indeed  this  is  how  old  systems  operated),  current  and  future 
processors  require  near-instantaneous  changes  in  threshold  levels  to  maintain  a  constant 
PFA  to  prevent  overloading  automatic  detection  and  tracking  systems  (Skolnik,  1980). 
One  such  processor  is  the  GO  CFAR  or  Greatest-Of,  Constant  False  Alarm  Rate  device. 
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B.  GO  CFAR  PROCESSOR 

In  the  GO  CFAR  processor,  the  background  noise  or  clutter  in  the  vicinity  of  a 
target  cell  is  sampled,  then  the  threshold  is  adjusted  to  provide  a  constant  false  alarm 
rate.  Figure  2  depicts  a  how  a  typical  radar  might  view  a  target  in  a  particular  region  in 
space. 


For  sake  of  simplicity,  assume  that  a  target  is  located  completely  within  a  single  range 
cell,  while  the  remaining  cells  contain  only  thermal  noise.  Once  received,  the  signal  is 
frequently  manipulated  to  create  a  quadrature  channel  which  is  90  out  of  phase  with  the 
in-phase  signal.  The  purpose  of  the  quadrature  channel  is  to  retain  the  phase  information 
and  eliminate  the  effects  of  blind  phases  (Skolnik,  1980). 

A  schematic  diagram  of  the  GO  CFAR  processor  is  shown  in  Figure  3.  The  input 
I  /  Q  thermal  noise  samples  are  considered  to  be  normally  distributed  N(0.,1)  and 
detected  as  either  x=  +  Q 2  or  x  =  a\I\+b\Q\  in  the  test  and  reference  cells.  Both 
reference  cell  neighborhoods  contain  n  cells  which  are  used  to  determine  the  noise  power 
levels  y„  y2.  The  detector  threshold  voltage  Vt  is  obtained  by  choosing  the  greater  of  y, 
and  y2,  normalizing  by  the  number  of  reference  cells  n,  and  multiplying  by  the  threshold 
multiplier  T.  A  target  is  declared  if  the  amplitude  of  the  test  cell  is  greater  than  Vr 
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Range  Gates - ^ 

X 


N(0,1)  _ _ _ 

Figure  3.  GO  CFAR  Processor  ( After  Skolnik,  1980) 

The  threshold  multiplier  (T)  is  normally  set  a  priori  to  achieve  a  desired  false 
alarm  rate  (Pace,  1994).  For  any  fixed  T,  increasing  the  number  of  cells  sampled  (#i) 
corresponds  with  a  decrease  in  the  false  alarm  rate.  The  probability  of  detection,  on  the 
other  hand,  is  based  upon  the  SNR  of  the  received  signal  and  is  normally  outside  the 
realm  of  the  designer's  control. 

For  mathematical  convenience,  CFAR  processors  are  frequently  analyzed  using 
the  envelope  detector  xe  =  Jl2  +  Q1  to  detect  the  In-phase  (I)  and  Quadrature  (Q) 
channels  of  the  baseband  signal.  To  reduce  hardware  complexities  created  by  the  use  of 

squares  and  square  roots  in  the  envelope  detector,  numerous  envelope  approximations 
have  been  developed  in  the  form  x  =  a\I\  +b\Q\  where  a  and  b  are  simple  multiplying 
coefficients  (Pace,  1994).  These  simplifications  result  in  a  slight  loss  in  sensitivity,  a 

difference  in  PFA  performance  (Pace,  1994)  and  consequently,  a  decrease  in  detection 
probability.  This  study  analyzes  the  detector  performance  of  the  GO  CFAR  processor 
using  the  envelope  detector  and  several  envelope  detector  approximations.  Since 
numerical  integration  of  the  PD  expressions  require  considerable  time,  Monte  Carlo 
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simulations  arc  used  to  obtain  the  results.  Closed  form  expressions  for  the  PD  are  then 
found  using  curve  fitting  techniques. 


C.  NOISE  AND  PROBABILITY 

Radar  signal  analysis  requires  an  initial  understanding  of  noise  which  is  a  random 
phenomenon. 


1.  Probability  Density  Functions 


Many  sources  of  noise  can  be  modeled  by  a  gaussian  or  normal  density  function, 
also  commonly  referred  to  as  the  bell  curve  as  illustrated  below  in  Figure  4, 


Figure  4.  Gaussian,  Normal  or  Bell  Curve 


and  given  by 


(x— T|)2 

f[x)  =  -\=e  2°2  (3.1) 

where  rpmean  and  o  =  standard  deviation,  N(r|,a).  The  x-axis  represents  signal 
amplitude.,  the  y-axis  the  probability  of  a  signal  associated  with  that  amplitude  Using 


this  model,  noise  is  comprised  primarily  of  low  amplitude  noise  centered  around  zero, 
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whereas  the  probability  of  higher  amplitude  noise  drops  off  exponentially  with 
increasing  amplitude. 

2.  Cumulative  Distribution  Functions 

The  cumulative  distribution  function  (cdf)  is  the  probability  that  the  value  x  is 
less  than  some  specified  value  and  is  given  mathematically  as  the  integral  of  the  density 
function 


Fiy)  =  T  ftx)dx. 

J  —oo 


(3.2) 


By  convention,  cdfs  are  represented  by  capital  letters,  while  pdfs  are  represented  by 
lower-case  letters.  The  cdf  of  the  density  function  in  Figure  4  appears  in  Figure  5. 


y 


Figure  5.  Typical  Cumulative  Distribution  Function. 
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3.  Summation  of  Density  Functions 

Given  two  different  functions  of  similar,  but  independent  random  variables  the 
pdf  of  the  output  is  equal  to  the  convolution  of  the  density  functions  (Papoulis,  1984): 

fziz)  =  f~  fx(z -y)fyiy)dy  (3.3) 

J  — oo 

Note  that  if  fx(x)  =  0  for  x  <  0  and  fy{ 0)  =  0  for  y  <  0,  then  fz{ 0)  =  0  for  z  <  0  and 

fz(z)=\zJAz-y)fy(y)dy  z>  0.  (3.4) 

4.  Maximum  of  Two  Density  Function 

Given  two  functions  of  independent  random  variables  fx(x)  and  fy(y),  cdf  of  the 
maximum  of  the  two  functions  is  given  as  (Papoulis,  1984): 

FAz)  =  FAz)Fy{z)  (3.5) 

and  the  corresponding  pdf  is  given  as: 

fz(z)  =  fx(z)Fy(z)  +fy(z)Fx(z).  (3.6) 

If  the  pdfs  of  fx(x)  and  fy(y)  are  equivalent,  then  the  pdf  simplifies  to 

fz(z)  =  2 fx{z)Fx{z)  =  2 \fy{z)Fy{z) .  (3.7) 

5.  Linear  transformation  of  density  functions 

For  the  linear  transformation  y =ax  where  x  is  a  random  variable  with  pdf  fx(x), 
the  pdf  of  y,  fy(y)  is  (Papoulis,  1984): 

1^4)-  <3-8» 
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6.  Absolute  value  of  a  density  function 

For  the  absolute  value  of  a  function  of  a  random  variable  y=|x|,  the  resulting  pdf 
is  given  as 

fy(y)  =/x(y)  +fx(-yl  0.9) 

The  relationships  discussed  above  are  needed  to  determine  the  density  function  for  the 
detector  approximation  being  examined. 

D.  ENVELOPE  DETECTOR  AND  APPROXIMATIONS 

The  input  to  the  GO  CFAR  processor  is  composed  of  a  combination  of  in-phase 

and  quadrature  signals  which  are  90°  out  of  phase  and  of  the  form  /  =  A  cos  (<j>)  +  x  and 
Q  =  A  sin  ((|>)  +  v  where  A  is  the  signal  amplitude,  <(>  is  the  signal  phase  and  x,  v  are 

Gaussian  random  variables  with  zero  mean  and  variance  of  one  N(0,1).  The  associated 
pdfs  of  the  absolute  value  of  these  two  signals  and  associated  curves,  with  <J>  =  71/8  are 
shown  below: 

(x-A  cos  f<t>))^  (-x-A  cos  (<t>))^ 

f\!\{x,A,§)  =  - - - - -j= - - -  X>0  (3.10) 


0  1^34 


Figure  6.  PDF  of  |I(  with  ({) — tt/8 
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(jC“v4sin(<l>))2  (-x-A  sin  (6)) 2 

f\Q\{x,A,§)  =  - - - — - - -  x  >  o  (3.11) 


Figure  7.  PDF  of  |Q|  with  (J)=p/8 


Using  Equation  (3.8)  to  factor  in  the  a  and  b  coefficients  for  the  envelope  approximation, 
we  obtain  the  following  pdfs: 


-^cos(<» 


fa\I\(x,§)  = 


|)2  [-i-.4cos(6)]' 

+e~  2 

\ci\j2K 


(3.12) 


[j-A  sin  (<J>)^ 2  sin  (it))]  2 

h  \Q\M-)  =  £=^Z^i  1  *>°  <313> 

As  previously  discussed,  the  envelope  detector  is  given  by  =  Jf  +  Q 2  while 
the  approximation  formula  is  x  =  a\I\  +  b\Q\.  Seven  sets  of  multiplying  coefficients  for  a 
and  b  which  provide  reasonable  approximations  to  the  envelope  detector  are  given  in 

Table  1  below. 


10 


Approximation 

a 

b 

i 

i 

1 

2 

1 

1/2 

3 

i 

1/4 

4 

i 

3/8 

5 

31/32 

3/8 

6 

0.948 

0.393 

7 

0.96043 

0.39782 

Table  1. 

Multiplying  Coefficients  a  and  b  for  Seven  Envelope  Approximations  (After  Filip,  1976) 

E.  PROBABILITY  OF  FALSE  ALARM  VS  THRESHOLD  MULTIPLIER 


Using  the  results  of  previous  sections,  the  PFA  of  the  GO  CFAR  processor  is 
shown.  The  GO  CFAR  processor  is  reproduced  in  the  figure  below. 


Input 
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-► 
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mit 
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x  =Jf+Q"  x  =a-  |lj  +b-  |q! 


Greater 

of 


y2 


Comparator 

— * — 
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N(0, 1 ) 


Figure  8.  GO  CFAR  (After  Skolnik,  1980) 


Since  we  re  interested  in  the  PFA,  we  start  with  the  assumption  that  no  target  is 
present  in  the  test  cell  and  ask  ourselves:  What  is  the  probability  that  the  signal  Tz/n  is 
higher  than  the  test  cell  output  xn  The  answer  to  this  question  gives  us  the  probability  of 
false  alarm.  We  start  with  the  pdf  of  the  input  signal,  step  through  the  GO  CFAR  model 
calculating  corresponding  pdfs  at  each  junction,  and  end  up  with  the  PFA  at  the  output  of 
the  GO  CFAR. 
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With  no  target  present  (A=0),  the  pdfs  of  |I|  and  |Q|  are  greatly  simplified.  In  fact, 
they  are  the  same. 

f\,\{x)=f\Q\  =  H[e^)  <3'4> 

For  the  envelope  approximation  x  =  u|/|  +  6|£?|,  we  use  our  previous  results  to 
calculate  the  pdf  of  the  input  as  the  convolution  of  the  pdfs  of  a|I|  and  b\Q\.  The  resultant 

pdf  becomes 


where  erf[\)  is  the  error  function  given  by 

er j{x)  =  -j=  Jq  e~~dy.  (3.16) 

Assuming  the  input  contains  only  noise  with  no  target  present,  Equation  (3.16)  also 
represents  the  pdf  of  a  test  cell  output. 

The  summation  of  input  signals  which  result  in  the  outputs  yl  and  y2  are  a 
function  of  the  number  of  reference  cells  n.  Using  the  results  of  Equation  (3.4).  we  note 
that  /’/(yl)  and  X,(y2)  are  both  n-fold  convolutions  of  /v(x)  for  which  there  is  no 
closed-form  solution. 

The  pdf  of  z  after  the  Greater-Of  block  can  be  determined  using  Equation  (3.7): 

f2(z)  =  2fyl  C z)Fyi  (z)  =  2 fy2(z)Fy2(z)  (3.17) 

where  Fw( z)  is  the  cumulative  distribution  function  offyl(z).  The  PFA  at  the  comparator 
output  is  then 
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PFA(T)  =  Jo"/z(z){  fnJAX)dk}dz. 


(3.18) 


Substituting  in  Equation  3.18  and  integrating  Fyl  we  get 

PFA(T)  =  2  J”/Vi  (3.19) 

where  fx(X)  is  given  Equation  (3.15).  Although  this  equation  has  no  closed  form  solution, 
it  can  be  calculated  numerically  Also,  closed  form  expressions  for  the  PFA  are  given  in 
Pace.  These  results  are  used  to  extract  the  corresponding  threshold  multipliers  used  in  the 
Monte  Carlo  simulations  for  the  detection  performance. 

F.  PROBABILITY  OF  DETECTION  VS.  SNR 

The  determination  of  the  probability  of  detection  follows  a  similar  line  of 
analysis  as  in  the  previous  section.  We  now  assume  however,  a  target  is  present  in  the 
test  cell.  In  addition,  the  threshold  multiplier  T  is  no  longer  a  variable,  but  a  constant 
which  corresponds  with  a  particular  PFA.  This  time,  we  ask  the  question:  What  is  the 
probability  that  the  signal  from  the  test  cell  containing  the  target  is  higher  than  the  signal 
Tz/n  appearing  at  the  comparator?  The  answer  to  this  question  gives  us  the  probability  of 
detection. 

Our  first  step  is  to  calculate  the  pdf  of  the  input  signal  x  using  the  envelope 
approximation.  Using  Equation  (3.4)  once  again,  we  can  calculate  the  pdf  of  the 
envelope  approximation  as  the  convolution  of  the  pdfs  of  a|I|  and  />|Q|  with  a  target 
present  (Equations  (3.12)  and  (3.13)  respectively): 


(3.20) 


Substituting  in  the  following  relation  of  SNR  to  signal  amplitude  A  (Wilson,  1982), 


(3.21) 


and  performing  the  integration we  obtain  the  following  closed-form  solution 


with  the  following  substitutions: 

A-  £ _ 

(3.23) 

2\a\\b\j2n(a2+b2) 

B  =  J2SNR  ( a  cos(<|>)  +  b  sin(<j>)) 

(3.24) 

C  =  JlSNR  (a  cos(<|))  -  b  sin(4>)) 

(3.25) 

D  =  a2+b 2 

(3.26) 

E  =j2SNRb2  cos® 

(3.27) 

F=  j2SNRahm(ty) 

(3.28) 
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G  =  JlSNRab  cos(<|>) 


(3.29) 


H=  J2SNR  ab  sin((|>) 

(3.30) 

l— ax 

(3.31) 

J  -  bx 

(3.32) 

The  PD  equation  is  similar  in  form  to  the  PFA  equation,  Equation  (3. 19),  but  T is 
now  constant  and  derived  for  a  specific  PFA. 


PD(SNR,ty)  =  J"  2/vl(z)(^/vl(^)4}[J”/n/,(>.,W,(|))ca]rfz  (3.33) 


Since  the  phase  of  the  target  signal  is  also  random,  the  PD  must  be  averaged  over  the 
phase  as 

PD(SNR)  =  \^^2fyi(z)\\^(%)^\j~dm,SNRX)^dzd^.  (3.34) 

Once  again,  there  is  no  closed  form  expression  for  the  PD  equation,  but  it  can  be 
evaluated  using  numerical  integration  or  Monte  Carlo  simulations.  The  latter  approach  is 
the  subject  of  the  following  chapter.  This  is  necessary  in  order  to  develop  the  closed  form 
expressions  to  approximate  this  detection  performance. 


15 


III.  MONTE  CARLO  SIMULATIONS  USING  MATLAB 


A.  MATLAB  PROGRAM 

Monte  Carlo  simulations  are  conducted  using  the  MATLAB  code  located  in  the 
appendix.  The  goal  is  to  obtain  the  PD  vs.  SNR  curves  associated  with  a  PFA  value  of 
10-4  for  each  envelope  approximation.  The  PFA  chosen  is  a  value  generally  reported  in 
the  literature  for  CFAR  comparisons.  For  both  the  envelope  approximation  and  envelope 
detector,  the  following  number  of  reference  cells  n  =1,2,4,8,16,32  are  used.  Threshold 
multipliers  (7)  for  each  particular  envelope  approximation,  the  number  of  reference  cells 
(«)  and  the  desired  PFA,  are  obtained  by  interpolating  the  PFA  vs.  the  threshold 
multiplier  data  from  a  previous  study  (Pace,  1994).  The  tables  below  summarize  the 
threshold  multipliers  used  for  each  Monte  Carlo  simulation. 


Table  2.  Threshold  Values  for  Envelope  Approximation  a=l,  b=l 


Table  3.  Threshold  Values  for  Envelope  Approximation  a=l,  b=l/2 
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I  n 

T  for  PFA=10"t 

1 

14.114 

2 

6.834 

4 

4.96 

8 

4.288 

16 

4.017 

|  32 

3.908 

Table  8.  Threshold  Values  for  Envelope  Approximation  a=0.96043,  b=0.39782 


n 

T  for  PFA=10'4 

i 

11.962 

2 

5.751 

4 

4.251 

8 

3.729 

16 

3.5215 

32 

3.4395 

Table  9.  Threshold  Values  for  Envelope  Detector  JF+Q* 


B.  PROGRAMMING  TECHNIQUES 

To  simulate  noise-only  inputs  to  the  reference  cells  ,  the  program  creates  matrices 
of  random  numbers  N(0,1)  with  dimensions  (n  x  N),  where  n  =  number  of  reference  cells 
per  side  and  N  =  number  of  simulations  desired.  The  absolute  values  of  these  matrices 
are  summed  across  the  n  variable  to  create  a  vector  of  dimension  (lxN),  then  multiplied 
by  the  appropriate  coefficients  a  or  b  and  added  together.  This  creates  two  vectors  of 
dimension  (1  x  N)  which  represent  the  summation  of  reference  cells  on  the  left  and  right 
side  of  the  test  cell.  These  two  vectors  are  compared,  and  for  each  N,  the  greater  of  the 
two  cells  is  placed  into  a  new  vector  z  of  dimension  (1  x  N).  To  create  the  comparison 
voltage  Vt  we  divide  by  the  number  of  reference  cells  n,  and  multiply  by  the  threshold 
multiplier  T.  The  test  cell  signal,  Vtest,  is  created  in  a  similar  manner,  but  this  time  we 
add  a  signal  with  amplitude  A  which  corresponds  to  the  SNR  to  be  tested. 

The  next  step  compares  the  two  vectors  Vt  and  Vtest  and  calculates  the 
percentage  of  times  that  Vtest  >  Vt.  For  each  SNR  tested,  the  process  is  repeated  for 
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0«jK7t/2  since  the  phase  of  noise  is  equally  distributed.  The  result  is  a  probability  of 
detection  for  a  single  SNR.  This  must  be  repeated  for  the  entire  range  of  SNRs  to  be 
tested. 

After  a  number  of  trial  runs,  it  was  determined  that  400,000  simulations  per  SNR 
with  an  interval  of  0.5  dB  provided  reasonably  smooth  curves  for  our  analysis.  We  also 
took  advantage  of  the  fact  that  at  some  point,  when  the  SNR  becomes  large  enough,  the 
PD  reaches  1.0.  Thus,  the  algorithm  breaks  the  loop  after  five  successive  PDs  of  1.0  to 
remove  unnecessary  computational  loops.  The  PD  plots  are  given  below. 
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Probability  of  Detection  Probability  of  Detection 


Probability  of  Detection 
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IV.  APPROXIMATION  FORMULAS  FOR  PD 

This  chapter  outlines  the  derivation  of  closed-form.,  curve  fit  expressions  for  the 
various  PD  vs.  SNR  curves. 

A.  METHODS  OF  CURVE  FITTING 

Numerous  methods  were  attempted  to  curve  fit  the  Monte  Carlo  results  and 
included  the  use  of  tanh,  third-order  exponentials,  and  ^//functions.  The  most  promising 
involved  the  use  of  the  erf  function. 

B.  ERF  FUNCTION  APPROXIMATION 

The  erf  or  error  function,  previously  given  in  Equation  (3.16),  is  repeated  here 
along  with  a  plot  of  the  function  for  convenience 


(4.1) 


-3  -2  -1  0  1  2  3 

X 


Figure  17:  Plot  of  Erf  Function 
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The  erf  function,  due  to  its  similar  appearance  to  the  PD  curves  in  the  previous  section,  is 
a  natural  candidate  as  a  curve  fit  solution.  By  rewriting  the  function  in  the  form 


f{x)  =  0-5(er^; 


x— cl 
c2 


+  1 


(4.2) 


we  can  control  both  the  curvature  and  axes  with  the  proper  values  of  cl  and  c2.  Using  the 
PD  vs.  SNR  curve  for  n=l  as  an  example,  we  can  obtain  a  reasonably  good  fit  using 
c=22.5  and  d=  5.26. 


8 

& 


P 


1c 


£ 

Cu 


PD  for  a=l,  b=l 
0.5(erf((x-22.5)/5.26)+l ) 


Figure  18:  Comparison  Between  Actual  PD  Curve  and  Curve  Fit  Solution  Using  Erf  Function 


Despite  the  symmetric  nature  of  the  erf  function,  we  are  still  able  to  produce  a  reasonable 
curve  fit  to  the  asymmetric  solid  PD  curve  (the  line  we  are  attempting  to  curve  fit)  and, 
at  the  same  time,  obtain  a  relatively  low  residual  across  the  entire  SNR  axis.  The  residual 
is  defined  as 
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Residual(SNR)  =  \PD(SNR) -flSNR)\ 


(4.3) 


and  describes  the  error  or  difference  between  the  two  curves.  The  residual  plot  is  shown 
below. 


In  this  case,  the  maximum  residual  equals  0.0218.  That  is,  for  any  given  SNR,  the  erf 
approximation  yields  a  PD  result  within  0.0218  of  the  actual  value. 

To  improve  this  result,  the  curve-fit  solution  is  split  into  two  parts:  PDs  >  0.5  and 
PDs  <  0.5.  A  curve  fit  solution  is  obtained  for  the  upper  portion  of  the  curve  above  a  PD 
of  0.5,  and  a  second  solution  obtained  for  the  lower  curve.  The  formulas  are  given  as: 

PD  <  0.5  PD  =  0.5(erfi^£i)  +  \)  (4.4) 

PD  >  ».5  PD  =  0.5(er/(^ff-)+\)  (4.5) 
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where  the  coefficients  cl,  c2,  c3  are  determined  experimentally. 

In  this  case ,  the  maximum  error  is  reduced  to  0  0116.,  clearly  a  better  solution. 
Plots  of  the  improved  solution  and  corresponding  residual  curve  for  a=l,  b=l,  n=l,  are 
given  in  the  following  figures. 
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Figure  20:  Comparison  of  Monte  Carlo  Simulation  vs.  Curve  Fit  Solution 


Figure  21:  Residual  of  Approximation 
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TABLES  OF  COEFFICIENTS 


The  following  tables  provide  the  associated  coefficients  for  each  of  the  seven 
envelope  approximations  and  envelope  detector. 


n 

cl 

c2 

c3 

i 

5.88 

4.655 

2 

4.494 

3.726 

4 

12.964 

3.93 

3.229 

8 

11.546 

3.807 

2.928 

16 

10.8 

3.677 

2.832 

32 

10.384 

3.602 

2.866  | 

Table  10.  Coefficients  for  PD  Curve  PFA-IO-4,  a—  1 ,  b—  1 
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cl 

c2 

c3 

n 

24.267 

6.901 

6.426 

2 

17.433 

5.311 

5.591 

4 

14.165 

4.643 

5.136 

8 

12.468 

4.34 

5.008 

16 

11.573 

4.191 

5.041 

32 

11.132 

4.14 

5.002 

Table  13.  Coefficients  for  PD  Curve  PFA=10,  a-1,  b— 3/8 


Table  14.  Coefficients  for  PD  Curve  PFA=10'4,  a— 3 1/32,  b— 3/8 


Table  15.  Coefficients  for  PD  Curve  PFA-10  ,  a-0.948,  b-0.393 


Table  16:  Coefficients  for  PD  Curve  PFA=10^,  a=0.96043,  b-0.39782 
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Table  17.  Coefficients  for  PD  Curve  PFA=10~4,  Envelope  Detector 


V.  PD  APPROXIMATION  COEFFICIENTS 

A  curve  fit  of  the  coefficients  as  a  function  of  the  number  of  reference  cells  n, 
allows  system  designers  to  extract  the  appropriate  coefficients  and  obtain  approximate 
PD  vs.  SNR  curves  for  any  number  of  reference  cells  n.  The  complete  set  of  coefficient 
plots  are  given  below. 


Figure  22:  Coefficient  Cl  as  a  Function  of  n 
For  Envelope  Approximation  a=l,  b=l 


Figure  23:  Coefficients  C2  (Solid  Line)  and  C3  (Dashed  Line)  as  a  Function  of  n 
For  Envelope  Approximation  a=l,  b=l 
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Coefficients  (c2  and  c3) 
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Figure  24:  Coefficient  Cl  as  a  Function  of  n 
For  Envelope  Approximation  a=l,  b=l/2 
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Figure  25.  Coefficients  C2  (Solid  Line)  and  C3  (Dashed  Line)  as  a  Function  of  n 
For  Envelope  Approximation:  a=l,  b=l/2 
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efficients  (c2  and  c3)  I  Coefficient  (cl) 


Coefficients (c2  and  c3)  |  |  Coefficient  (cl) 


Figure  28:  Coefficient  Cl  as  a  Function  of  n 
For  Envelope  Approximation  a=l,  b=3/8 


Figure  29.  Coefficients  C2  (St  m  and  C3  (Dashed  Line)  as  a  Function  of  n 
For  Envelop. .  . iximation:  a=l,  b=3/8 
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Coefficient  (cl) 


Figure  30:  Coefficient  Cl  as  a  Function  of  n 
For  Envelope  Approximation  a=31/32,  b=3/8 


Figure  3 1 .  Coefficients  C2  (Solid  Line)  and  C3  (Dashed  Line)  as  a  Function  of  n 
For  Envelope  Approximation:  a=31/32,  b— 3/8 
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Coefficients  (c2  and  c3)  I  Coefficient  (cl) 


Figure  36:  Coefficient  Cl  as  a  Function  of  n 
For  Envelope  Detector 


Figure  37.  Coefficients  C2  (Solid  Line)  and  C3  (Dashed  Line)  as  a  Function  of  n 

For  Envelope  Detector 
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VI.  CONCLUSIONS 


Monte  Carlo  simulations  yield  accurate  PD  vs.  SNR  curves  for  envelope 
detection  approximation  GO  CFAR  processors  These  curves  can  be  closely 
approximated  with  closed  form  expressions  using  the  erf  function  in  the  form: 

pd  <  0.5  PD  =  0.5(art^4  +  1)  (6.1) 

PD  >  0.5  PD  =  0.5(er/(^f-)  +  1 )  (6.2) 

where  the  coefficients  cl.,  c2  and  c3  are  chosen  for  a  specific  envelope  approximation 
and  number  of  reference  cells  n.  These  PD  approximations  compare  favorably  with  the 
results  of  Monte  Carlo  simulations;  the  maximum  residual  in  all  cases  is  less  than  0.028. 

A  curve  fit  of  the  coefficients  cl,  c2  and  c3  as  a  function  of  the  number  of 
reference  cells  provides  a  means  of  quickly  extracting  the  coefficients  for  any  number  of 
reference  cells  1<«<32  for  a  PFA=10'4.  In  conjunction  with  the  equations  above,  these 
coefficients  produce  accurate  PD  curves  as  a  function  of  SNR. 
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APPENDIX  -  MATLAB  PROGRAMS 


Programs  to  Run  Monte  Carlo  Simulation 

runmont.m  Batch  file  to  run  multiple  simulations 

setvar.m  Program  to  set  variable  for  each  loop 

montesim.m  Main  program  to  run  simulations  of  envelope  approximations 

monteenv.m  Main  program  to  run  simulation  of  envelope  detector 


Programs  to  run  Plots  of  Data 
plotmont.m  Plots  data  from  *.dat  files 
getdata.m  Gets  data  to  plot  graphs 
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RUNMONT.M 


%  FILE  NAME:  RUNMONT.M 

%  Batch  file  to  run  all  Monte  Carlo  Simulations  for  PFA=le-4 
clear 

%  Input  parameters 
%  a:  coefficient  a 

%  b:  coefficient  b 

%  n:  number  of  cells  per  side 

%  PFA:  Required  PFA 

%  T:  Threshold 

%  N:  number  of  monte  carlo  simulations  per  inner  loop 
%  MAXSNR:  Maximum  SNR  in  dB 

% - 

%  Set  variables  using  data  from  pfa_a*.dat  at  PFA=le-4 
%  SNR/PD  data  saved  to  file  MCl__4.dat 
set__var;  %  Clear  all  variables  and  reset  to  clear  memory 
a-1;  b=l; 

nTABLE=[  1.000  2.000  4.000  8.000  16.00  32.000]; 
TTABLE=[  12.287  5.952  4.423  3.891  3.680  3.5975]; 
montesim; 

save  MCl_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_b*.dat  at  PFA=le-4 
%  SNR/PD  data  saved  to  file  MC2_4.dat 
set_var;  %  Clear  all  variables  and  reset  to  clear  memory 
a=l;  b=0.5; 

nTABLE=[  1.000  2.000  4.000  8.0000  16.00  32.000]; 
TTABLE=[  13.430  6.515  4.773  4.1525  3.903  3.8040]; 
montesim; 

save  MC2_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_c*.dat  at  PFA-le-4 
%  SNR/PD  data  saved  to  file  MC3_4.dat 
set_var;  %  Clear  all  variables  and  reset  to  clear  memory 
a=l;  b=0.25; 

nTABLE=[  1.000  2.000  4.000  8.0000  16.00  32.000]; 
TTABLE^t  16.685  7.915  5.530  4.6750  4.335  4.197]; 
montesim; 

save  MC3_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_d*.dat  at  PFA=le-4 
%  SNR/PD  data  saved  to  file  MC4_4.dat 
set_var;  %  Clear  all  variables  and  reset  to  clear  memory 
a=l;  b=3/8; 

nTABLE==[  1.000  2.000  4.000  8.0000  16.00  32.000]; 
TTABLE=[14.530  7.018  5.064  4.3610  4.076  3.9636]; 
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montesim; 

save  MC4_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_e*.dat  at  PFA=le-4 
%  SNR/PD  (L;a  saved  to  file  MC5_4.dat 
set_var,  %  Clear  all  variables  and  reset  to  clear  memory 
a=31/32;  b=3/8; 

nTABLE=[  1.000  2.000  4.000  8.0000  16.00  32.000]; 
TTABLE=[14.399  6.960  5.032  4.3380  4.058  3.9465]; 
montesim; 

save  MC5_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_P.dat  at  PFA=le-4 
%  SNR/PD  data  saved  to  file  MC6_4.dat 
set_var;  %  Clear  all  variables  and  reset  to  clear  memory 
a=0.948;  b=0.393; 

nTABLE=[  1.000  2.000  4.000  8.0000  16.00  32.00]; 
TTABLE=[14.1 12  6.835  4.960  4.2870  4.017  3.907]; 
montesim; 

save  MC6_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_g*.dat  at  PFA=le-4 
%  SNR/PD  data  saved  to  file  MC7_4.dat 
set_var;  %  Clear  all  variables  and  reset  to  clear  memory 
a=0.96043;  b=0.39782; 

nTABLE=[  1.000  2.000  4.000  8.0000  16.00  32.00]; 
TTABLE=[14.1 14  6.834  4.960  4.2880  4.017  3.908]; 
montesim; 

save  MC7_4.dat  mctemp  -ascii 

% - 

%  Set  variables  using  data  from  pfa_P.dat  at  PFA=  1  e-4 

%  SNR/PD  data  saved  to  file  MCr_4.dat 

set_var;  %  Clear  all  variables  and  reset  to  clear  memory 

nTABLE=[  1.000  2.000  4.000  8.0000  16.000  32.000]; 

TTABLE=[1 1.962  5.751  4.251  3.7290  3.5215  3.4395]; 

monteenv; 

save  MCr_4.dat  mctemp  -ascii 
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SETJVAR.M 


%  FILE  NAME:  SET.VAR.M 

%  File  to  set  variables 

clear  global 
clear  all 

%  Input  parameters 
%  a:  coefficient  a 

%  b:  coefficient  b 

%  n:  number  of  cells  per  side 

%  PFA:  Required  PFA 

%  T:  Threshold 

%  Nmax:  Max  number  of  monte  carlo  simulations  per  angle  phi 
%  MAXSNR:  Maximum  SNR  in  dB 

angle=pi/2;%  Set  max  angle  to  sample  (0  to  pi/2) 

phi=0:angle/40:angle;%  Set  phi  to  vector  of  angles  from  0  to  pi/2,  40  steps 
PFA=le-4; 

Nmax=10000; 

MAXSNR=40; 
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MONTESIM.M 


%  FILE  NAME:  MONTESIM.M 
%  Generic  Monte  Carlo  Simulation 
% 

%  Input  parameters 
%  a:  coefficient  a 

%  b:  coefficient  b 

%  n:  number  of  cells  per  side 

%  PFA:  Required  PFA 

%  T:  Threshold 

%  Nmax:  Max  number  of  monte  carlo  simulations  per  angle  phi 
%  N:  Size  of  random  matrix  (nxN)=Nmax 
%  MAXSNR:  Maximum  SNR  in  dB 

%  nTABLE:  Values  of  n=number  of  cells  per  side 

%  TTABLE:  Values  of  T=corresponding  threshold  setting  for  given  n  and  PFA 
%  angle:  Max  angle  to  sample  (0  to  pi/2) 

%  phi:  vector  of  angles  from  0  to  pi/2,  40  steps 


for  num=  1 :  length(  nT  ABLE) ; 

skip=0;  %  Sets  flag  to  zero  (Used  to  remove  unnecessary  loops) 
count=0;  %  Sets  count  to  zero  (Used  to  remove  unnecessary  loops) 

n=nT  ABLE(  num) ; 

T-TTABLE(num); 

N=round(Nmax/n);  %  Size  of  random  matrix  (nxN)=Nmax 

SNRdB=0:0.5:MAXSNR;%  Set  Signal-to-Noise  ratio  in  dBs  (increments  of  0.5  dB) 
PD=zeros(  1  ,length(SNRdB)); 


% - 

for  i=  1 : length( SNRdB ) ;  %  Performs  loop  over  all  SNRs 
if  skip==l 
break 
end 


for  j=l:length(phi),  %  Performs  loop  over  all  angles  phi:  0  to  pi/4  (20  steps) 

SNR=1 0A(SNRdB(i)/l  0);%  Converts  SNR  from  dBs  to  linear  scale  [1,1] 


A=sqrt(2*SNR); 

Inl=randn(n,N); 

In2=randn(n,N); 

Qnl=randn(n,N); 

Qn2=randn(n,N); 

if  n=l 


%  Converts  SNR  value  to  required  A  [1,1] 
%  1  on  left  side,  noise  only  [nxN] 

%  I  on  right  side, noise  only  [nxN] 

%  Q  on  left  side,  noise  only  [nxN] 

%  Q  on  left  side,  noise  only  [nxN] 


y  1  =a*abs(In  1  )+b*abs(Qn  1 ); 
y2=a*abs(In2)+b*abs(Qn2); 
else 

yl=a*sum(abs(Inl))+b*sum(abs(Qnl));  %  Sum  of  left  side  noise  inputs  [lxN] 
y2=a* sum( abs( In2 ) )+b * sum( ab s( Qn2 ) ) ;  %  Sum  of  right  side  noise  inputs  [lxN] 
end 
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z=max(yl,y2);%  GO  output  [lxN] 

Vt=T*z/n;  %  Comparison  voltage  [lxN] 

% - 

%  Test  voltage  in  test  cell  with  signal  [1,N] 

%  z=a*|I|+b*|Q|  or  z=a*|Acos(phi)+x|+b*|Asin(phi)+y| 

Vtest=a*abs(  A*cos(phiO'))+randn(  1  ,N))+b*abs(  A*sin(phi(j))+randn(  1  ,N)); 

% - 

%  Calculate  percentage 
%  Add  up  all  percentages  for  each  value  of  phi 

PD(i)=PD(i)+length(find((Vtest-Vt)>0))/(N*length(phi));  %  [lx#  of  SNRs] 
end 

%  deletes  unnecessary  loops  if  PD=1  for  greater  SNRs 
if  PD(i)>— 1 
count=count+l ; 
if  count>5 
skip=l; 
end 
end 
end 

% - Break  branches  to  here 

for  k=  1 :  length!  SNRdB) , 
if  k>=i 
PD(k)=l; 
end 
end 

%  Store  SNR/PD  data  to  file  named  mctemp  (temporarily) 

mctemp((2*num)-l:2*num,:)=[SNRdB;PD]; 

end 
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MONTEENV.M 


%  FILE  NAME:  MONTEENV.M 
%  Generic  Monte  Carlo  Simulation  for  envelope  detector 
% 

%  Input  parameters 
%  n:  number  of  cells  per  side 

%  PFA:  Required  PFA 

%  T:  Threshold 

%  Nmax:  Max  number  of  monte  carlo  simulations  per  angle  phi 
%  N:  Size  of  random  noise  matrix  (nxN)=Nmax 
%  MAXSNR:  Maximum  SNR  in  dB 

%  nTABLE:  Values  of  n=number  of  cells  per  side 

%  TTABLE:  Values  of  T=corresponding  threshold  setting  for  given  n  and  PFA 
%  angle:  Max  angle  to  sample  (0  to  pi/2) 

%  phi:  vector  of  angles  from  0  to  pi/2,  40  steps 


for  num=  1 :  length!  nTABLE); 
skip=0;  %  Sets  flag  to  zero 
count=0;%  Sets  count  to  zero 
n=nTABLE(num); 

T=TTABLE(num); 

N=round(Nmax/n);  %  Set  size  of  random  noise  matrix  (nxN)=Nmax 


SNRdB=0:0.5:MAXSNR;%  Set  Signal-to-Noise  ratio  in  dBs  (increments  of  0.5  dB) 
PD=zeros(  1  ,length(SNRdB)); 

% - 

for  i=l:length(SNRdB),  %  Performs  loop  over  all  SNRs 
if  skip=l 
break 
end 


forj=l:length(phi),  % 

SNR=  1 0A(  SNRdB(  i)/l  0); 
A=sqrt(2*SNR); 
Inl=randn(n,N);  % 

In2=randn(n,N);  % 

Qnl=randn(n,N);  % 

Qn2=randn(n,N);  % 


Performs  loop  over  all  angles  phi:  0  to  pi/4  (20  steps) 
%  Converts  SNR  from  dBs  to  linear  scale  [1,1] 
%  Converts  SNR  value  to  required  A  [1,1] 

I  on  left  side,  noise  only  [nxN] 

1  on  right  side,noise  only  [nxN] 

Q  on  left  side,  noise  only  [nxN] 

Q  on  left  side,  noise  only  [nxN] 


ifn— 1 


yl=sqrt((Inl  A2)+(Qnl.A2)); 
y2=sqrt((In2.A2)+(Qn2.A2)); 
else 

yl=sum(sqrt((Inl  A2)+(Qnl.A2)));  %  Sum  of  left  side  noise  inputs  [lxN] 

y2=sum(sqrt((In2.A2)+(Qn2.A2)));  %  Sum  of  right  side  noise  inputs  [lxN] 

end 

z=max(yl,y2);  %  GO  output  [lxN] 

Vt=T*z/n;  %  Comparison  voltage  [  1  xN] 

% - 

%  Test  voltage  in  test  cell  with  signal  [1,N] 
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%  z=sqrt[(IA2)+(QA2)  or  z=sqrt[(Acos(phi)+x)A2+(Asin(phi)+y)A2] 
Vtest=sqrt((A*cos(phi(j))+randn(  1  ,N))  A2+(A*sin(phi(j))+randn(l  ,N))  A2); 

%  Calculate  percentage 
%  Add  up  all  percentages  for  each  value  of  phi 

PD(i)=PD(i)+length(find((Vtest-Vt)>0))/(N*length(phi));  %  [lx#  of  SNRs] 
end 

%  deletes  unnecessary  loops  if  PD=1  for  greater  SNRs 
if  PD(i)>=l 
count=count+l ; 
if  count>5 
skip=l; 
end 
end 
end 

% —  Break  branches  to  here 

for  k=l  :length(SNRdB), 
if  k>=i 
PD(k)=l; 
end 
end 

%  Store  SNR/PD  data  to  file  named  mctemp  (temporarily) 

mctemp((2*num)-l:2*num,:)=[SNRdB;PD]; 

end 
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PLTMONT.M 


%  FILE  NAME:  PLTMONT.M 

%  Batch  file  to  plot  PD  vs.  SNR  curves  from  Monte  Carlo  Simulations 

%  Plot  PD  vs  SNR  curve  from  mcl_4.dat 

loadmcl__4.dat 

save  mc.dat  mcl_4  -ascii 

varjitle-PD  vs.  SNR  a=l,b=l,PFA=le-4,  n=  1,2, 4, 8, 16,32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc2_4.dat 

load  mc2_4.dat 

save  mc.dat  mc2_4  -ascii 

var_title='PD  vs.  SNR  a=l,b=l/2,PFA=le-4,  0=1,2,4,8,16,32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc3_4.dat 

load  mc3_4.dat 

save  mc.dat  mc3_4  -ascii 

var_title='PD  vs.  SNR  a=l,b=l/4,PFA=le-4,  n=  1,2, 4, 8, 16, 32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc4_4.dat 

load  mc4_4.dat 

save  mc.dat  mc4_4  -ascii 

var_title- PD  vs.  SNR  a=l,b=3/8,PFA=le-4,  n=  1,2, 4, 8, 16, 32’; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc5_4.dat 

load  mc5_4.dat 

save  mc.dat  mc5_4  -ascii 

var_title='PD  vs.  SNR  a=31/32,b=3/8,PFA=le-4,  0=1,2,4,8,16,32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc6_4.dat 

load  mc6_4.dat 

save  mc.dat  mc6_4  -ascii 

var_title=’PD  vs.  SNR  a=0.948,b=0.393,PFA=le-4,  0=1,2,4,8,16,32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mc7_4.dat 

load  mc7_4.dat 

save  mc.dat  mc7_4  -ascii 

varjitle=’PD  vs.  SNR  a=0.96043,b=0.39782,PFA=le-4,  n=  1,2, 4, 8, 16, 32'; 
get_data 

%  Plot  PD  vs  SNR  curve  from  mcr_4.dat 

load  mcr_4.dat 

save  mc.dat  mcr_4  -ascii 

var_title='PD  vs.  SNR  for  envelope  detector, PFA=le-4,  n=  1,2, 4, 8, 16, 32'; 
get_data 
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GET_DATA.M 


%  FILE  NAME:  GET__DATA.M 

%  Batch  file  to  plot  PD  vs.  SNR  curves  from  Monte  Carlo  Simulations 

%  Get  data 
elf 

load  mc.dat 

num  _plot=length(mc( : ,  1  ))/2 ; 
hold  on 
grid  on 

fori=l:nurn_piot; 

plot(mc((2*i)-l,:),mc((2*i),:)) 

end 

xlabelfSignal-to-Noise  Ratio  (dB)*) 
ylabel(’Probability  of  Detection’) 
title(var_title) 

%  print  -dwin 
pause 
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